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ABSTRACT 

The use of superconductors in the field windings of large 
synchronous rotating electrical machines virtually eliminates the 
need for magnetic iron core material and hence reduces the size and 
weight requirements for a given power level. This possible reduction 
in size and weight can be exploited in the propulsion of hydrofoil 
c y a ^* Present hydrofoil propulsion schemes all have serious salient 
disadvantages that have stimulated the search for a better propulsion 
scheme for hydrofoil craft. This thesis proposes the use of a novel 
dual armature superconducting induction motor. The motor would be 
encased in a water tight pod and placed on the end of a foil strut 
and directly drive a super-cavitating propeller. A mathematical 
model is developed and the size and parameters of a 21000 horsepower 
motor are determined. The resulting model is statically and 
dynamically analyzed. 
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I INTRODUCTION 

In conventional displacement type ship hull forms, the power avail- 
able within the hull can be readily converted to propulsive thrust by 
conventional direct drive propellers. However, in a hydrofoil type ship 
hull form, the mechanism of converting the available power within the hull 
to propulsive thrust is no longer a simple matter. When the hydrofoil 
vehicle 'flies' (becomes foil borne), the hull is lifted out of the water 
and conventional types of propulsion are no longer feasible. A simple 
direct drive straight rotating shaft with a propeller on the end will 
not suffice. 

An assortment of alternative hydrofoil propulsion schemes have 
been proposed or developed: 

1. The water jet system scoops up water at an inlet in a forward 
foil strut and ejects this water at high velocity out the stern, 
thus providing the propulsive thrust. This arrangement requires 
large amount of internal hull volume for water ducts and pump 
machinery and places large amounts of water weight inside the hull. 
Additionally, this method has a low efficiency at low velocities. 

2. Another scheme utilizes a mechanical gear train containing 
right angle gears to drive propellers mounted on the foil struts. 

This scheme suffers from a very high initial cost and a low mean 
time between failure, hence unreliability and high maintenance costs. 

3. Air propellers, turbo-props, and other aerodynamic propulsors 
all have los efficiencies in the speed ranges anticipated for 
hydrofoils in the near future. 
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4. Electric motors can be placed in pods on the end of the foil 
struts and drive propellers directly. This method eliminates the 
objections of the preceding schemes, but introduces a new 
objection of its own. The pod of the size required for a conventional 
motor would produce extremely high levels of hydrodynamic drag in 
the speed ranges of the hydrofoil vehicles. 

While all of the above methods of propulsion are feasible and can be 

utilized, there is much needed improvements before hydrofoil vehicles 

can compete economically with the displacement type vessels. Much 

research and development is being done to improve on each type of propul~ 

sion system. In particular, the motor in the submerged pod system has 

been vigorously pursued with emphasis on reducing the required size of the 

motor through the use of super conduction technology. 

The first attempt at applying super conducting materials to motors 

( 7 ) 

utilized Faraday's homopolar concept. It has the characteristics of 
low DC voltages with high current densities .^* ^ The high DC current 
densities of the homopolar devices complicates the power transmission and 
switching problems and requires some sort of liquid brushes in the motor. 
Several different configurations using the homopolar concept are 
discussed in reference (4). 

The recent trend in super conducting machines has been toward the 
use of AC synchronous devices. While AC machines lack the high starting 
torques which are characteristic of DC machines, they have a superior 
reliability in a marine environment and the power transmission and 
switching problems are much less than for DC machines. The AC super 
conducting machines most considered are of the type with a stationary three 
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phase armature (stator) and a revolving DC field (rotor) consisting of 
super-conducting winding in a Dewar bottle. 

A recent development in AC super-conducting machines is a family of 
machines which utilizes two armatures and a super-conducting field 
winding. The field winding and one armature* which may or may not be 
wound, rotate, and the other armature is fixed (stator). One motor config- 
uration utilizes a superconducting field which is free to rotate and is 
mechanically unconstrained except for its own inertia, it responds only 
to the applied electrical torques and its inertia. Both armature windings 
are wound in a three phase manner, one armature is fixed (stator) and 
receives the input electrical power, and the other, positioned between the 
stationary armature and the rotating field, develops the mechanical output 
torque in the same manner as a wound rotor induction motor does. This 
latter configuration is the topic of this thesis. It will be mathematically 
developed and applied to the propulsion of a hydrofoil vehicle. Its static 
and dynamic characteristics will be presented. The resulting size of the 
motor will demonstrate that it is compact and will significantly reduce the 
hydrodynamic drag of the pod to an acceptable level. 
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II PHYSICAL DESCRIPTION OF THE DUAL ARMATURE MOTOR 

The motor to be studied by this thesis utilizes a four pole 

rotating field that is mechanically unconstrained, except for its own 

moment of inertia. This field will be of the super conducting type and 

will be encased in a Dewar bottle and cooled by liquid helium. It will 

be excited by a DC current source. In general, the design of the field 

( 18 ) 

will be similar to the field winding developed at MIT. This field 

will be capable of producing conventionally unachievable magnetic flux 
densities due to the large current densities that are made possible by 
the super conducting material. The superconducting winding will hence 
forth be referred to as the field. 

Surrounding the field and attached to it is a damper shell. This 
shell will provide the start up torques and damn out dynamic oscillations. 
The damper shell and the field together will hence forth be referred to as 
the rotor. 

Surrounding the rotor will be a three phase, four pole winding. 

This winding will contain no magnetic core material and will contain 
only enough steel or other material to provide mechanical support and 
rigidity. It will be free to rotate relative to the rotor and to the 
motor frame. It will be attached to a shaft from which mechanical power 
may be extracted and supplied to the load, in this case a propeller. The 
electrical terminals of the winding will be brought outside the motor 
through slip rings. This winding will hence forth be referred to as 
the shell. 

The rotor and shell are surrounded by yet another three phase four 
pole winding that contains no magnetic core material. It is 



rigidly attached to the motor frame and contains only enough steel to 
provide mechanical support and rigidity. This winding will hence forth 
be referred to as the stator. 

Surrounding the stator is a shield made of a magnetic core material 
Its purpose is to con-tain the magnetic fluxes within the motor and 
provide a more efficient return path for the fluxes. If the shield were 
not present the rotating flux field would extend beyond the motor and 
possibly induce hazardous voltages in the surrounding areas which could 
present a hazard both to equipment being operated in the vicinity and 
to man. 

The stator is now surrounded by a wa.ter tight enclosure , called the 
pod, that will keep the motor dry and minimize the hydrodynamic drag. 

The magnetic shield could be a part of the pod body. 

In operation, three phase electrical power will be provided to 
the stator winding which will produce a revolving magnetic field. The 
rotor will loci: into synchronism with the stator field and in effect 
become a virtual magnetic core by producing a rotating magnetic field 
far greater than could the air cere stator by itself. The shell will 
act as a wound rotor induction armature under the influence of the 
rotating magnetic field and drive a super cavitating propeller directly 
as shown in figure (1). Figure ( 2 ) shows the geometric configuration 
of the field, damper, shell stator and shield. 
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(4) 3-**6i4 
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Figure ( 2 ) 
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III MATHEMATICAL MODEL OF THE MOTOR 

The motor is modeled electrically as shown schematically by figure 
(3)* Figure ( 3 ) shows the axes of the various windings. The rotor has 
a direct axis (in the direction of positive flux lines of the field) and 
a quadrature axis which leads the direct axis by ninety degrees. The 
direct axis is at an angle 0 with respect to the stator 'A* phase 
axis as shown in figure (3), 

The shell 'A* phase axis is at an angle ^ with respect to the stator 
•A' phase axis. 

The various mutual inductances are assumed to be sinusoids dependent 
upon the angular orientation of the axes involved. The mutual inductances 
are maximum when the axes are aligned, zero when they are ninety degrees 
apart, and negative (minimum) when they are one hundred-eighty degrees 
apart. Equation set (1) gives the flux linkages for the various circuits 
when the currents in all circuits are considered. Equation set (2) gives 
the voltages at the terminals of the various circuits when all currents in 
all circuits are accounted for and the circuit relation *= p\^ + R^ 1^ 
is utilized. The symbol p is a operator representing the time derivative 
d/dt. All other symbols utilized are idetified in the Glossary of Terms, 
Table (1). 

The rather formidable equation sets (1) and (2) can be significantly 
simplified by performing a Park's transformation on them. The Park's 
transformation equations are given here as equation set (3). When the 
transformation of equation set ( 3 ) is applied to equation sets (l) and 
(2), equation sets ( 5 ) and (6) result. 
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Table I Glossa.ry of Terms 

Notes where a dual assignment of usage is indicated, the context should 
make it clear to the reader the intended usage of the symbols. 

Symbols 
I Current 

L Self Inductance 

1 length 

M Mutual Inductance 

N Coil Turns/ angle as defined in figure ( 3 ) 

P Permeance 

p Differential Operator d/dt 

R Resistance/ Radius 
s Slip 

T Torque 

t time 

U Dummy Variable used in the Park's Transformation Equation 
V Voltage 

W Angular Velocity 

X Reactance/ Ratio of Reactances 

0 Angle Defined in Figure ( 3 ) 

X Flux 

tt 3.n4.... 

/ Angle Defined in Figure ( 3 ) 



Per-Unitized Flux 



Table I 



continued 



lk 



Subscripts 

a Phase a/ Stator Armature 

b Phase b/ Base Quantity 

c Phase c 

d Direct Axis 

f Field 

i Arbitrary Index 

j Defined in Text 

k Defined in Text 

1 Defined in Text 

m Defined in Text 

o Zero Sequence/ Synchronous Frequency/ Defined in Text 

p Per-Unitized Quantity 

r Rotor 



s 



Shell 
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CkCK 




Figure (3) 
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Equation set (4) can now be obtained directly from equation sets (5) 
and ( 6 ). 

A further simplification may be made to equation set (5) by 
per unitizing it in accordance with reference (3). First* equation set (5) 
is re-written as equation sets ( 7 ), ( 8 ) and ( 9 ). 

Then* the ordinary variables in equation sets (?), ( 8 ) and ( 9 ) are 
normalized in accordance with reference ( 3 ) by dividing each variable 
quantity by a corresponding base quantity. The results are shown as 
equation sets ( 10 ), ( 11 ) and ( 12 ). 

It is highly desirable that the resulting per-unitized mutual 
inductances exhibit the property of reciprocity* that is = M^. 

This imposes the following restrictions on the base quantity selection. 



2. v I = 2. y I 

2 ab ab 2 sb sb 



V I = V I 
rb rb fb fb 



Additionally* the following conditions on the base quantities are 
arbitrarily imposed to simplify the per-unitized equations. 

Vi I , 0 M W I v M W I (L-M ) W I , 

„ _ 3 as o sb _ 3 ar o rb af o fb _ a aa ' o ab 

x a ” 2 — vT — = — vT — " — vT — vT 

ab ab ab ab 

Thus, at this point, six constraints have been placed on the eight base 
quantities. This leaves only two base quantity parameters free to be 
chosen arbitrarily • Or.ce these two more conditions are specified the base 
quantities will be completely determined. Further consideration of base 
quantities is given in appendix (H). 

Equation sets (10), (11) and (12) may now be written as equations sets 



(13), (14) and (15) 
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Appendix (a) shows that X 



( s £ and X^ = X ^ so equation 



sets (13)» (14) and (15) can be re-written as equation sets (16), (17) 
aiid (18) respectively. 

If the base quantities utilized to obtain equation sets (10), (11) 
and (12) are applied to equation set (4), equation set (I9) results. 

Now,, if equation set (19) is solved for the time derivatives of 
the fluxes, equation set (20) results. 

Equation sets (16), (17) and (18) can be inverted so that fluxes 
are the dependent variables and currents are the independent variables. 
The resulting equations are shown as equation sets (21), (22) and (23); 

where, _ X s ( X s ” X r } 

X j " X ( X - X J 
d a v a s' 

( X - X ) 

' C* 7 



h m TT 

' ra 



X ) 
s 



X 1 - 



( X - X ) 
'a r 



( X - X ) 
y = ' B V ' 

x m rxr-Xn 



A o Tx“"x~y 

' r s' 

Equation sets (21), (22) , (23) and (20) can be combined to obtain 
equation set (24). 

The direct and quadrature axis voltages of the stator armature are 
V 



adp 



V COS ( 0 ) and V 
a v aqp 



- V a SIN ( 0 ) 



given by; 
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The dynamic mechanical equations of the shell and rotor in per 

unit form are: p 2 ( 9 ) - W T /z H 

o rp r 

and p 2 (0-N)= W (T -T. )/Z H 

r v ' o K sp lp 7/ s 

Appendix (b) provides the following expressions: 



T = 7 ^ ( '■£' '-P - 'i' ) 

rp W ' rd sq sd rq ' 



sp W ^ sd aq ad sq ' rp 

These equations along with the equations of equation set (24) form 
a complete set of differential equations that mathematically describe 
the motor's model. It is assumed in this thesis that all loadings and 
transient states will be symetrical and balanced. Thus, v iP and 

3LO 

so will be assumed to be zero. Additionally, ^ ^ will be considered 
a constant. This being the case, the given differential equations, 
reduce to a ninth order non-linear set of differential equations. 

These equations can be solved once all the constant parameters are 
known. This solution is rather straight forward for the steady state 
condition, although a computer minimises the computational effort. 

The transient state defies a hand solution. The digital computer 
can be used to simulate the transient state once a set of initial 
conditions has been determined. 

References ( 8 ), (9)» (10) and (12) provided valuable background 
material that aided in the development of the equation in this section. 



IV Determination of Machine Parameters 
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The size of the motor is determined in Appendix ( C ) and the circuit 
parameters are determined in Appendix ( D ). The results are summarized in 
Table ( 2 ) . 



Table 2 



Machine Size and Circuit Parameters 
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Inner Radius of Field 


4.0 inches 


Outer Radius of Field 


6.5 inches 


Average Radius of Damper 


?.0 inches 


Inner Radius of Shell 


7.5 inches 


Outer Radius of Shell 


10. 5 inches 


Inner Radius of Stator 


11.0 inches 


Outer Radius of Stator 


13*5 inches 


Inner Radius of Shield 


14.5 inches 


Outer Radius of Shield 


20.0 inches 


Electrical Length 


7.4 feet 


R ap 


0.20 


R 

sp 


0.588 


R rp 


26.87 


% 


0.0 (super-conducting winding) 


X. 

3 


7.799 


X k 


3.685 


X 1 


4.685 


X 

0 


10.08 


H 

s 


0,738 seconds 


H 

r 


0.0738 seconds 
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V Steady State Characteristics 

When the motor is running in steady state, the governing 
differential equations of equation set ( 24 ) can be reduced by setting 
all time derivatives of the fluxes to zero. Additionally, the torque 
on the rotor must be zero and this directly implies T = SP *» 0 . 

Slip, s, is defined defined in the conventional manner; that is, 

W - W 

s = — - where = P ( 8 ) “ W Q and W g “ p ( ) 

w 0 - p ( 4 ) 

in the steady state, then s = . As can be seen from 

o 

figure ( 3 ) e 0 = N + ^ so p( H) = W q - p ( ^ ) , so finally p (N) = s W q . 

Taking equation set ( 24 ) and setting the time derivatives of 

fluxes to zero, letting “ 'P =• 0 , and using the relations 

pN = s W q , p0 = W q and I_ = / R , equation set ( 25 ) results. 

The torque equation of the shell reduces to, 

X 

<p a «. \Jj 

sp W Q sd aq and by using the fourth equation of 

equation set ( 25 ), the expression for torque can be written as, 

R X 'P 0 
T « ( _o )2 

sp s v W J 

r o 

Solving for 'F yields " if ~ E " 

^ ^ sp v sp 

Assuming that slip, s, stator torque, T , and field current, I„, 

2 p X 

are the independent variables, then equation set ( 25 ) can be solved 
for the fluxes and voltages which result in equation set ( 26 ). 

The digital computer is a valuable tool in evaluating the expressions 
of equations set ( 26 ) numerically. Equation set ( 26 ) and the machine 



34 



R X 

0 - V , + + a P - (X.'P.-X.'P.) 

adp aq W Q v k sd j ad ' 



R X 

0 - V - ¥ . + - —ff — ( - X . ¥ ) 

aqp «id W© v j aq ' 



0 = X. ^ - X- ** + ^ , 

k ad 1 sd rd 



0 = s ¥ + - X. ¥ 

sd W k aq 
o 



°”' i 'sd -(‘ +X „ ) s k + X m' i ’f 



0 - I f + ( * -'*% > 



rd * f 



Equation Set ( 25 ) 



35 






axl 




W o 



aq X. X 
^ k o 



W - 0 
sq 



R X 

V£i S - P ' . ° . x 'i* 

^sd s W A k aq 

o 



*i» =0 

rq 






W o X f 



rd sd X 



W 

VX( t= up + _J? T 

f rd X X x f 

o m 



R X 

V = _ »j> + — 2P — £ ( y if 

v adp r aq W q ^ A j Y ad 



X. SP, ) 
k sd ' 



R X 

I - SE» . + -- a jy ■■ X. 'P 
aqp ad W j aq 



V, = /(v . ) 2 + (V ) J 
t V v adp' v aqp' 



Equation Set ( 26 ) 



parameters listed in Table ( 2 ) along with the load characteristics 
determined in appendix ( E ) were injected into an IBM model 370 
computer use the WATFIV program listed in appendix ( G ), 

The results of the computer analysis are shown in the results 
section as figure ( 4 ) through ( 1 5 )• 



VI Transient State Characteristics 
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Equation set ( 24 ), table ( 2 ) and the shell and rotor torque 
expressions of Appendix ( B ) provide the mathematical discription 
of the transient state. Given these equations and a complete set of 
initial conditions and driving functions, the system is deterministic. 

The initial conditions are determined in Appendix ( F). 

The transient state was simulated in four distinct evolutions: 

1 . Rotor start up with R = infinity 

sp 

2. Shell start up with rotor at full speed and at various values of 
R 

sp 

3. Crash back. 

b. Three phase fault at the stator armature terminals. 

The simulations were all conducted using an IBM model 370 digital 
computer and the IBM supplied program CSMP (Continuous Systems Modeling 
Program). The program listings are shown in Appendix ( G ). 

The rotor start up simulation assumed that the shell circuit is 
open and in effect not there. This was equivalent to starting a 
conventional single armature machine. The results of this start up 
ar shown in figure ( 16 ). 

The shell start up simulation was conducted by first having the 
rotor running in steady state then placing various values of finite 
resistance in the shell circuit. The results of this simulation are 
shown in figure ( 17 ). 

The crash back simulation first assumed a steady state full power 
condition for the motor, then the phase sequence of the stator terminal 
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voltage were reversed. The results of the crash back simulation are 
shown in figure ( 18 ). 

The three phase fault simulation again assumes a steady state full 

power condition for the rotor, then the three terminal leads of the stator 

armature are electrically shorted. In effect V = V , = V =0. 

aa ab ac 

The results of the three phase fault are shown in figure ( 19 )• 
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VII Results 



Figure ( 4 ) through figure ( 15 ) show the steady state analysis 
results. In these figures, s, 1^ and T are the independent variables. 
Figures ( 4 ), ( 5 ) and ( 6 ) show the terminal voltages. 

Figures ( 7 ), ( 8 ) and ( 9 ) show the terminal currents. 

Figures (10), (11) and (12) show the power factor. 

Figures (13), (14) and (15) show the rotor angle. 

Figure ( 16 ) through figure ( 19 ) show the angular velocities of the 
shell in the transient state. 

Figure ( 16 ) shows the rotor start up transient. 

Figure ( 1? ) shows the stator start up transient. 

Figure ( 18 ) shows the crash back transient. 

Figure (19) shows the three phase fault transient. 
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Figure ( 5 ) 



Armature Terminal Voltage pu 
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Figure ( 6 ) 



Armature Terminal Current pu 










Figure ( '/ ) 
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Figure ( 9 ) 
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Figure ( 12 ) 
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Figure ( 14 ) 
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Angular Velocity of Shell and Rotor 
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VIII DISCUSSION OF RESULTS AND CONCLUSIONS 

Figures (4) through (15) demonstrate that the dual armature 

motor can be characterized by a set of steady state condition curves. 

These curves can assist in the rational selection of an operating 

condition for the motor. The data supplied by these curves can be 

cut many different ways; for instance, I„ or T could be considered 

fp sp 

as dependent variables and or 1^ could be treated as independent 

variables. In cutting the data different ways additional data points 

may be desirable, in this case the WATFIV program in Appendix (G) could 

be used if the appropriate parameters were adjusted. 

Figures (16) through (19) provide the transient state results. 

Figure (16) shows that the rotor on start up under goes large angular 

velocity excursions with the attendant large acceleration forces. At 

one point the rotor actually runs backwards. Eventually the rotor 

settles down to its mechanical synchronous speed of 1800 RPM. 

Figure (17) shows that for the shell start up there is some critical 

parameter R between 10 pu and 20 pu where the shell starts up in a 

stable fashion at the fastest possible rate. For values of R less than 

sp 

this critical value, the rotor tumbles or pulls out of synchronism with the 

stator field and results in a loss of starting torque and the introduction 

of a pulsating torque. For R = O .588 pu (its minimum value) the motor 

sp 

rotor does not start at all but oscillates about the zero point. 

Figure (18) shows the crash back W curve. The results of this 
curve beyond 1.8 seconds are questionable. The computer simulation appears 
to have gone astray at this point, perhaps due to too large of a step 
size in the simulation integration routine or due to the proximity of W o 
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to zero. The rotor and shell both move at approximately the same velocity 
in this maneuver. 

Figure (19) shows the three phase fault at the stator armature 
terminals. Again the computer simulation appears to go astray at about 
3.7 seconds so results beyond that point are questionable. 

The preliminary results of this thesis indicate that the dual 
armature super conducting motor is indeed a candidate for hydrofoil 
propulsion and warrants further investigation. 
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IX RECOMMENDATIONS 

Additional areas that could be investigated using the equations 
developed in this thesis nay be: 

1. Determine steady state conditions as a function of R . The 

sp 

present investigation considered R = O .588 pu in the steady state 
condition. 

2. Write a computer program that will accept machine dimensions as 
input and return machine parameters for use in the steady state 

or transient analysis. 

3. Instead of holding ^ constant as done in this thesis in the 
transient state investigation, hold I^ constant. 

4. Investigate current densities in the various windings during 
the transient periods to see if they exceed tolerable levels. Devise 
a scheme for controlling excess currents during transients. 

5 . Determine the static stability limit of rotor angle. 

6. Refine load characteristics. 

7 . Vary machine dimensions to optimize the design. Criteria for 
optimising could also be developed as a side issue. 

8. Introdue hydrofoil dynamics into the model. 

9. Investigate speed control by varying W q the electrical fre- 
quency of the prime mover. 

Most of these subject areas would apply at an under graduate project 



level 
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Appendix A Shell Constraints 



Equation sets ( 13 )» ( 14 ) and ( 15 ) contain parameters that are 
not independent. The following development assumes that flux 
leakages between adjacent circuits is negligible. 

Assuming that the stator, shell and field axes are aligned so as 
to eliminate the sinusoidal dependency of the mutual inductances, the 
following relations can be written. 

< L s - H ss > ‘ 2 N s 2 < P + 6 P ss > 

" N r 2 < P + 6 ?r > 



M = N N ( P + 6 P ) 
as a s v as ' 



H 



M 

M 



ar 

af 



M 

M 



sr 

sf 

rf 



= N N ( 
a r v 



- "a N f < 



= N N ( 
s r v 



tl s H f ( 



H r N f < 



P + 
P + 
P + 
P + 
P + 



6 P 



ar 



6 P 



af 



6 P 



sr 



6 P 



sf. 



6 P „ 
rf 



) 

) 

) 

) 

) 



The factor appearing in the first equation is characteristic of 
three phase windings. 

Now, if the various leakage permeances may be considered negligible 
and are set to aero, the above equations may be used to derive the 
following relations. 
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2 . m Cl _m ) = 2. m * 2 m 
2 af ' a n aa ' 2 sf 2 as 



4 H (L-K ) =■ ^ M * 4 M 
2 ar v a as. 7 2sr 2 as 



M . L = M - M 
af r rf ar 



These equations may now be per-unitized as follows s 



0 M . W I , (L - M ) W I, , M W I, 0 M W I . 

2 . af o ab v s ss' o sb _ 2 sf o sb 2 as o ab 



fb 



sb 



fb 



sb 



- M WI.(L-M)WI. 0 M W I . - M W I , 

2 asr o ab v s ss 7 o sb 2 sr o sb 2 as o ab 



rb 



sb 



rb 



sb 



M q v; I L W I M . W I M tf rt I 

ax o fb r o r d r f o f b ar o rb 



V 



V 



ab rb rb ab 

Comparing these relations with those in equation sets ( 10 ) 
( 13 ) it can be seen that the above relations reduce to; 



X s " X sf, 

X - X . 
s sr and 



X = X _ 
r rf . 



and 
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Appendix B 



Derivation of Torque Expressions 



The torque on any winding in per unit is given by 



T ip ■ < * iq T idp * ^id 



Applying this to the stator and shell windings and substituting the 
expressions for current from equation Sets (21) and (22) and canceling 
terms where possible yields: 



T 1 S3 

ap W 



X X T 
o k 



( ‘F . - ‘f 1 ¥ . ) and 

' ad sq aq sd 



T <= — — ( vp - qJ ¥ ) 

sp VJ o ^ T sd rq J ' sq T rd ' 



ap 



Also equilibrium of forces requires T + T + T =0, 

^ u ap sp rp 

Solving the above equations for the torque on the rotor gives 
X 

T = — ( ‘F p _ <p p ). 
rp W v sq rd T sd T rq ' 



Appendix G Kotor Sizing 
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References (11) and (14) contain equa.tions for the power rating of 
an air core, three phase synchronous motor. They can be simplified and 
adapted to the present application by assuming the following} 

a. number of pole pairs = 2 

b. winding angle of field « 120° 

c. winding angle of shell *= 60° 

d. winding angle of stator <= 60° 

e. all winding electrical lengths= 1 

f. the following parameters are defined as 



X = R ai/ R ao * y = R fi/ R fo » and z " • The radius 



si' so 



dimensions axe shown in figure (C-i). 

The resulting power rating expression is, 



p - ‘V 6 10 "' 1 “o J a J f c 2 V - jo< r > 



where O, - | ( -ln.t ( i ~r sL ) ( f) , 

sh 



- a - X a sin(y) = - 7 gos 2 (y) , 



COS(y) = power factor , 

h V °° 75 



J * G s , \o x4 
J f C 2 I 1 ' ^ ( H fo ’ 



c s . i ( x 4 1„ X + + U .. -x u l\ !k f ). 



sh 



and 
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Figure ( G - 1 ) 
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The current densities were chosen based on the values used in 

reference (11) and are; J 13 2.50 x 10 amps/meter , 

8 2 

and 53 1.25 x 10 amps/meter . 



Synchronous electric frequency was chosen to be W Q = 377 radians/sec. 
Power factor was chosen to be PF = 0.8 . 

The power was chosen based on the requirements specified in reference 
(4) and is P = 21,000 horespower per motor. 

The radii selected are; 



R . =11.0 inches 
ai 



R^ = 7.0 inches 



R^ = 4.0 inches 

R . =7*5 inches 
si ^ 



R„ =6.5 inches 
fo 

R = 10.5 inches 
so 

R = 13® 5 inches 
a,o 

R gh = 14.5 inches 



Finally, if the above assumptions are made and the power rating 



expression is solved for 1, the result is 1 = 7.4 feet 
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Appendix D Circuit and Mechanical Parametex-s of Motor 

References ( 11 ) and ( 14 ) contain equations for determining the 
inductances (both self and mutual) for an air core, three phase, 
synchronous motor. These equations are infinite series summations; 
however, the lead term, n=l, dominates and all other terms are small. 
These equations can be simplified and adapted to the present application 
by assuming, in addition to the assumptions of appendix (c), that; 

a. the leading term of the infinite series dominates, 

b. the shield is made of magnetic iron and 

c. the theory of superposition holds . 

The resulting inductance expressions are listed in Table ( D - 1 ). 

The dimensions are the same as discussed in appendix (c) and shown in 

figure ( C-l ). If these dimensions are now used to solve the expressions 

in table ( D-l ), table ( D-2 ) results. 

If the expressions of table ( D - 2 ) are per unitized using the 

base quantities of appendix ( H ), table ( D-3 ) results. 

The definitions of X^, X^, X^, X m and X q given on page 25 can now be 

evaluated and are 5 X. ® 7*799 

X k = 3.685 

X x » 4.685 

X <= 10.08 
o 

Reference (14) provides an equation for R^ and is 

WIN 2 

r £_ 

r 4 a ft. 

s t 

where O g = o t^, t^ being the damper shell thickness and a being the 



Table D-l Inductances 



1 8 % 1 n 

^ (1-X 2 ) 2 L 



ln(x) + 



(1-X 4 ) (1-X 4 ) 2 

— — + — 

I- 8 



fao ,rj 
R sh -1 



18 u o ^ 1 □ 

’ (1 _ Z V L_ 



(i-z c y 



(1-Z 4 ) (1-Z 4 ) 2 , B so v 

ln(Z) + + ( ) 

4 8 R sh 



9 



IT u N 2 1 r "* 
L = o_x__ 

r 16 



1 + ( 



~ )4 1 
R , ~J 
sh 



L - 



27 u 0 1 



f 2 n 8 (l-Y 2 ) 2 



Y u ln(Y) + 



(1-Y 4 ) (1-Y 4 ) 2 , R f0 7 



8 



( — ) 

R sh 



M *» -|L 
aa a 



(1+Z 2 ) R, 



u N N 1 



so 



as ^3 ° a s (i_x 2 ) R 



( — r 



ao 



(1-X 4 ) R 



ao 



-ln(X) + 



— -( — r 

4 R sh - 



3 u N N 1 R 
J o a r 



oar / 

M ” 5 ( 

2TT (I"**) R 



ar 



i-fj: 



ln(X) + 



( !ao ,rj 



ao 



sh 



M 



af 2rr 3 oaf (1 _ x 2) 



E- 



(1-X 4 ) R 
( 

R 



ao 



sh 



4 



Table D-l Inductances ( cont. ) 
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H = 
sr 



3 u N N 1 
o s r 

2 TT (1-Z 2 ) 



R 



) 2 Qln(Z) + 



SO 



(l-z 4 ) 




M 



9 \/T 



U N N„ 1 



(1+y 2 ) r 



sf 2 „3 O 8 f - (1 . z2) ■ 



(^ 0)2 



so 



(i-z 4 ) R 



so 



-ln(Z) + — ( -sa ) 



R. J 

sh 



M - 



3vr 



* 16 it ^ W(^ 2 >(^)£ + (^7] 
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Table D~2 Inductances 

L = 5.104 x 10" ? N 2 1 

cl cl 

L « 3.532 x 10“ 7 N 2 1 
s s 

L = 2.601 x !0~ 7 N 2 1 
r r 

L f = 2.101 x 10" 7 N f 2 1 



Table D-3 Per Uniti&ed Inductances 
X = 9.089 

cl 

X s = 19.23 
x = 56.62 

r 



X = 112.4 
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conductivity of copper. Assuming a - 3*6 x 10^ m ^ OS/ ^meter and 
tj = 0.11 indies, R^ = 4.42 x 10 ^ N ^ 1. Per unitizing using the 
base quantities of appendix ( H ) gives, R = 26.87 pu. 

In general, the resistance of a wound three phase winding is 
given by 



R i = 



4 N. 2 1. 
1 1 



(R . 2 - R. . 2 ) e . K . o. 

v 01 ii ' wei si i 



assuming the conductivity is O - 6 x 10^ mhos/meter and the space factor 

K g = 0.27 for both the shell and the stator armature, the per unitized 

resistances are; R = 0.200 pu and R = 0.588 . 

ap s p 

The resistance of the field winding is of course zero because 
it is made of a superconductor. 

An initial estimate of the moments of inertia for the rotor and 

for the shell with shaft, propeller and entrained water included were ; 

J = 120 slug ft 2 and J = 12 slug ft 2 
s ^ r ° 

For H defined in the conventional manner, that is, 

2 

J. V 

u — -t- h 

H i 2 



- 0.738 sec and H « 0.0738 sec. 

Note: the computer listings of appendix (g) show H to be different 

from the above values. This is because W was included so the value in 

o 

the listings is 377 times smaller than the actual assumed value. 



Appendix E Load Characteristics 
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The dynamics of the motor and propeller system are 

considered to be uncoupled from the hydrofoil dynamics for the 

time span of the simulations since the time constants of the two 

systems are quite different. Thus the load torque characteristics 

are a function of W alone. 

s 

For the shell start up, the propeller characteristics are 



assumed to be 



bp ‘ 5 O' A) 2 - 



For the crash astern and the three phase fault, the load 
characteristics are assumed as shown below. 




/ 
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Appendix F Determination of Initial Conditions 

An operating point was selected for the full power steady state 
from figures (4) through (15) as s=0.00l5, 1^=1. 0 pu, and T g p=1.0 pu. 
Then using the equations of equation set (26) a complete set of initial 
condition were determined and are listed in the computer listings in 
appendix (G). 

For rotor start up, it was assumed that T =0.0, W =0.0, W =0.0. 

sp r s 8 

and 0=0.0 . The remaining initial conditions were calculated using 
Equation Set (26) and are listed in the computer listings in appendix 
(G). 

Part of the rotor start up computer simulation routine was to 
print out the final values of all the time variables at the end of 
the simulation period. These final values were then utilized as the 
initial conditions for the shell start up. The values are listed in 
the computer listings in appendix (g). 



Appendix G 



Computer Ftogram listings 
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Appendix H Base Quantity Selections 

As mentioned on page ( 21 ), two remaining base parameters are free 
to be selected. The bases selected are; 

1) = 21*000 horsepower and 

2) = leO pu when = 1,25 x 10® amps/meter^. 

At this point eight conditions have been specified on the eight 

base quantities so the base system is fully constrained and all base 
quantities can be calculated. They are; 

c 2,696 x i 0 5 aunps/N a 

^st = 4.4-79 x 1 0 5 amps/N s 
ly-^ “ 1,388 x 10® amps /N r 
= 2,217 x 10® amps/Np 



- 19»36 N a volts 
v _ II.65 „ .. 

V sb " N s volts 

v- _ 5*640 „ _ 

V rb “ N r voits 

vo it s 
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